Three-body recombination at large scattering lengths in an ultracold atomic gas 
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We study three-body recombination in an optically trapped ultracold gas of cesium atoms with 
precise magnetic control of the s-wave scattering length a. At large positive values of a, we measure 
the dependence of the rate coefficient on a and confirm the theoretically predicted scaling pro- 
portional to a 4 . Evidence of recombination heating indicates the formation of very weakly bound 
molecules in the last bound energy level. 



Three-body recombination, the process of two particles 
forming a compound while interacting with a third parti- 
cle, is of general importance in many-body systems. Un- 
derstanding three-body processes has for a long time been 
a challenge in many areas of physics. In dilute ultracold 
atomic gases, where unique experimental access to inter- 
particle interactions is available, three-body recombina- 
tion leads to the formation of diatomic molecules. The 
possibility to trap and cool such molecules and thus to re- 
alize ultracold molecular ensembles 1] or atom-molecule 
mixtures Q holds great prospects for future research in 
the field of quantum gases. 

The main experimental challenge for measuring three- 
body recombination is to distinguish three-body losses 
from two-body losses. In magnetic traps, this is in gen- 
eral difficult, as dipolar relaxation is present as an inher- 
ent two-body loss mechanism. In the special case of mag- 
netically trapped 87 Rb, experiments could nevertheless 
provide reliable data because of the anomalously weak 
two-body decay in this species 0, ■ As a more general 
approach, optical traps allow to store atoms in the low- 
est internal state where two-body loss is fully suppressed 
for energetic reasons 0. Three-body recombination in 
bosonic systems with large s-wave scattering length a 
has been explored in experiments utilizing Feshbach res- 
onances in 23 Na and 85 Rb Q. The results demon- 
strate the dramatic enhancement of loss processes near 
such resonances. 

In the case of large positive a, theoretical studies 
Hi EJ fiot HH predict formation of molecules in a weakly 
bound s level and deduce scaling of the three-body re- 
combination rate with a to the fourth power. Recombina- 
tion is discussed in terms of the event rate v IOC = a rec n 3 
per unit volume and time, where n denotes the parti- 
cle density. In the low-energy limit, the universal scal- 
ing law for the characteristic parameter a roc is given as 
aw = Cha 4 /m y|, where m is the mass of the atom, and 
C a dimensionless factor. Predictions for C give values 
between and ~70 @, 0j El EI , with oscillatory behav- 
ior between C = and C max = 67.9 expected for strong 
variations of a on a Feshbach resonance El- Previously 
available experimental data on three-body recombination 
for various species roughly agree with the general trend of 
a universal a 4 scaling, with evidence for deviations from 



this scaling law near a Feshbach resonance [6j. 

In this work, we employ an optical trap to measure 
three-body recombination in an ultracold thermal gas 
of 133 Cs in the 6 2 S 1 / 2 ,L 1 = 3,tof = 3 absolute ground 
state. In this state, the scattering length varies strongly 
with an applied magnetic field B through a combination 
of broad and narrow Feshbach resonances at low field 
strengths [l^ . Il3j | . The values of a as a function of B 
have been calculated to high precision El- This tun- 
ability allows us to conduct experiments in a range of 
precisely controlled large positive values of a, where the 
condition of recombination into a weakly bound s level 
is well fulfilled. Detailed measurements of the time evo- 
lution of both atom number and temperature allow us to 
quantify both three-body loss and heating. 

The three-body loss rate coefficient L3 is related to a TCC 
via -L3 = nia Tec , where ni denotes the number of atoms 
lost from the trap per recombination event. Therefore, 
the central relationship for comparing measurements to 
theory is 



L 3 = niC—a 
m 



(1) 



In the recombination process the molecular binding en- 
ergy e is set free as kinetic energy. The molecule and the 
third atom receive e/3 and 2e/3, respectively. Since usu- 
ally e is large compared to the trapping potential depth, 
both are expelled from the trap, setting ni = 3. How- 
ever, the binding energy of the weakly bound last energy 
level of the dimer is given by e = ft / (ma 2 ) El> and at 
very large scattering length values e may be below the 
trap depth and the third atom cannot escape. If the po- 
tential of the atom trap does not confine the molecule, 
the dimer is lost and ni = 2. If, however, the molecule is 
trapped and stays within the atom cloud, it may quickly 
quench its high vibrational excitation in an inelastic col- 
lision with a fourth trapped atom. The large amount of 
energy released in this situation expels the collision part- 
ners, making ni = 3. In cither case, the kinetic energy of 
the remaining atom is distributed in the ensemble, giving 
rise to recombination heating. 

Our measurements are conducted in a nearly isotropic 
optical dipole trap formed by two crossed 100 W CO2- 
laser beams focused to a waist of 600 ^m 16j. A vertical 
magnetic field gradient of 31 G/cm cancels gravity for 
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the 6 2 Si/2,-F = 3,tof = 3 absolute ground state of Cs. 
The optical trapping forces are so weak that no other 
rap substate can be held against the gravitational force, 
thus perfect spin polarization of the sample is enforced 
by the trap itself. We have measured the trap frequencies 
by exciting small oscillations along the principal axes, the 
geometric average results to u/(2tt) — 14.5±0.5 Hz. The 
trap depth is fcs x 12 /iK; fee denotes the Boltzmann 
constant. Trap heating is extremely weak, from long- 
time measurements we can give an upper bound for the 
heating rate of 0.2 nK/s. The background gas limited 
1/e lifetime is ~180 s. 

After loading the dipole trap 0]> the ensemble cools 
down by plain evaporation within the first 10 s to 
3 x 10 6 atoms at T = 1 fiK, 1/12 of the trap depth. 
To avoid further evaporative loss and cooling in our mea- 
surements, we cool the sample well below this limit. This 
is achieved by forced radio frequency evaporative cooling 
in the magnetic levitation field. We apply a radio fre- 
quency resonant with the Zccman splitting between two 
niF substates. In the vertical levitation gradient this 
acts as height selective removal from the trap, effecting 
ID evaporative cooling. Evaporation proceeds at a bias 
field of B = 75 G to tune the scattering length to 1200 do 
(do = 0.053 nm denotes Bohr's radius), which at the low 
initial peak density of 3x 10 11 cm -3 keeps the elastic scat- 
tering rate at ^200 s _1 , high enough for efficient cooling. 
A simple 7 s two-segment linear sweep of the radio fre- 
quency provides 1.3 x 10 6 atoms at 450 nK (peak phase- 
space density J}w4x 10~ 3 ). Some measurements have 
been taken at 200 nK starting temperature; to reach this, 
we appended another linear 5 s rf ramp, leaving 4 x 10 5 
atoms in the trap (D as 10~ 2 ). Temperature and atom 
number are measured by taking an absorption image with 
typically 50 ms of expansion in the levitation field after 
switching off the CO2 lasers. For the expansion, the bias 
field is in all measurements set to 75 G to ensure constant 
conditions. 

In a first set of measurements, we have recorded the 
atom number and temperature of the trapped sample af- 
ter a holding time of 10 s at fixed magnetic fields varied 
between 10 G and 100 G [13 ■ Figure [T] shows the data 
in comparison to the calculated scattering length [lij . 
which includes Feshbach resonances involving s- and d- 
wave channels. Additional resonances involving closed 
g-wave channels are marked by arrows. The Feshbach 
resonances have been observed in previous experiments 
through light-induced two-body losses 0, 0] and de- 
scribed theoretically [I2I lisj. In our measurements, all 
resonances appear as prominent loss and heating fea- 
tures. Over the full range, both atom number and tem- 
perature show a clear correlation with the magnitude of 
the scattering length. 

To obtain more quantitative results on recombination 
loss and heating, we have studied the time evolution of 
the atom number N and temperature T at various mag- 
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FIG. 1: Remaining fraction of atoms N/No (upper curve) and 
temperature T (middle) after 10 s in the C02-laser trap in a 
range of magnetic bias fields. Initially iVo = 1.3 x 10 6 atoms 
are prepared at 450 nK. For comparison, the lower plot shows 
the calculated scattering length a with additional Feshbach 
resonances indicated by arrows (see text). 
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FIG. 2: Atom number and temperature as a function of time 
in the trap at a bias field of 56 G. The solid lines show a 
numerical fit to the data. The corresponding three-body loss 
coefficient is L 3 = 7(2) x 10" 25 cm 6 /s. 



netic fields. Figure shows a typical measurement at 
B = 56 G. Integrating the local three-body loss rate L^n 3 
over the sample gives N/N — — L% (« 2 ). Expressing the 
average density in terms of the directly accessible quan- 
tities N and T, we replace L 3 (n 2 ) by "/N 2 /T 3 , where 
2 
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L3(mo; 2 /27r/cB) 3 /v / 27. With an additional loss term 



—aN for background gas collisions, we obtain a differen- 
tial equation for the atom number, 



dN N 3 
— = -aN- 1 ^. 



(2) 



Atom loss through three-body recombination leads to 
"anti-evaporation" heating. With its n 3 dependence, 
three-body recombination predominantly happens in the 
region of highest density in the trap center. The mean 
potential energy of an atom undergoing a three- 
body process can be calculated by integrating over a 
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FIG. 3: Measured values of L3 as a function of the bias field 
B. Filled symbols represent measurements taken at an initial 
temperature of 200 nK, open symbols at 450 nK. The solid line 
shows the calculation according to Eq. Q with n;C = 225. 
The scale on the right hand side displays the corresponding 
scattering length values. The inset shows the region around 
the 48 G resonance with a model fit (see text). 



thermal distribution with weights proportional to n 3 , 
which yields ^ksT. As the ensemble average is ffcaT, 
for each lost atom an excess energy of lk B T remains in 
the sample. To account for recombination heating, we 
introduce an additional constant energy fc^T/j per lost 
atom. Relating the total heating energy &b(T + Th) 
to the average energy 3k B T of a trapped particle yields 
f/T = (N/N) k B (T + T h )/(3k B T), and we obtain an 
expression for the temperature evolution, 



dT 
~dt 



1- 



N 2 (T + T h ) 



(3) 



T 3 3 

Equations and © form a set of two coupled nonlin- 
ear differential equations for TV and T. Using an iterative 
approach |20j , we do a least-squares fit of the numerical 
solutions of both equations to the experimental data. For 
the data in Fig. [2J the fit is plotted as solid lines. Fit pa- 
rameters are a, 7 and Th- Error estimates for the fit 
parameters are based on varying one parameter up and 
down, resp., to the points where x 2 increases by one, 
while optimizing the other parameters ■ 

The data obtained for L3 as a function of magnetic 
field is plotted in Fig. |3J The inset shows an expanded 
view of data taken on the 48 G Feshbach resonance. The 
error bars represent the statistical errors derived from the 
fit as described above, which typically range around 20%. 
The systematic uncertainty is larger: Due to the scaling 
of L 3 with N 2 u> 6 , small errors in these values amount to 
large differences in the result obtained for L3. Taking 
the uncertainty in the trap frequency and an estimated 
30% in atom number, the systematic error can be up to 
a factor of two. 
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FIG. 4: Recombination heat Th as a function of the bias field 
B. Filled symbols represent measurements taken at an initial 
temperature of 200 nK, open symbols at 450 nK. The inset 
shows the region around the 48 G resonance. The shaded area 
shows the expected range based on the calculated binding 
energy of the last bound state of the CS2 molecule (see text). 



In order to compare our data to theory, we do a least- 
squares fit of the expected scaling law [Eq. dl to our 
measurements, using the calculated a(B) Single fit 
parameter is niC. The resulting curve is drawn as a solid 
line in Fig. The data clearly confirm the universal a 4 
scaling. In our experimental range, where a varies by 
less than a factor of two, we do not find any indication 
of resonant behavior in C j^. The result from the fit 
is niC — 225, corresponding to C — 75 for nj = 3 or 
C = 112 for ni = 2. In comparison to theory, both values 
are somewhat above the upper limit of C max ~ 70, but 
in good agreement within our systematic error limit. 

The inset in Fig. 03 shows the data on the 48 G Fesh- 
bach resonance. In order to model the shape of the loss 
feature, we have to take into account the width of the 
sample in the vertical B field gradient and the unitarity 
limitation at the finite temperature. The theory curve 
is convoluted with a gaussian of full 1/e width 260 mG, 
which is given by the vertical extension of the n 3 dis- 
tribution of the trapped sample. A cutoff in a accounts 
for the unitarity limitation; the cutoff value of 1800 a® 
results from a least-squares fit to the data. Our _B-field 
scale is slightly adjusted for an optimum match with the- 
ory by introducing a 40 mG shift, which is well within 
our calibration accuracy. 

Via the condition k « 1, the inverse wavenumber 
k^ 1 = \AkBTm/ {ttT% 2 )\ _1 / 2 characterizes the crossover 
to the zero-energy limit. At 450 nK, k^ 1 = 1500 ao; 
at 200 nK, fc" 1 = 2250 a Q . Correspondingly, the 450 nK 
data on the Feshbach resonance show a strong effect from 
the unitarity limitation, while the data taken at 200 nK 
starting conditions follow the expected scaling law up to 
the maximum a sa 1500 ao. 

Our results on recombination heating are displayed in 
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Fig. 0] A trend towards smaller Th at higher magnetic 
fields, corresponding to higher scattering lengths a, is vis- 
ible. At large values of a close to the Feshbach resonance, 
the values are correspondingly very small. We compare 
the data to the scenario that the recombination heating 
arises just from the binding energy e of the last bound 
state of the dimer. As the atom involved in the three- 
body recombination process receives 2e/3, the recombi- 
nation heat per lost atom results to ksTh = 2e/(3n;). 
Therefore, we assume ksTh to be between 2e/9 (ni = 3) 
and e/3 fa = 2). This range is shown as shaded area in 
Fig. E] 24] . This simple model shows reasonable agree- 
ment with the experimental data. Our heating measure- 
ments thus provide evidence of the formation of very 
weakly bound ultracold molecules. 

From the apparent effect of recombination heating, we 
can draw an important conclusion for the efficiency of 
evaporative cooling towards quantum degeneracy in gases 
with large positive scattering length. If temperatures 
are very low (fc^T -C e) and the atom that carries 2e/3 
cannot immediately escape from the trap, recombination 
heating has a detrimental effect much worse than mere 
three-body loss. For efficient evaporation it is thus im- 
perative to avoid the hydrodynamic collision regime and 
to implement a three-dimensional evaporation scheme. 
In our cooling experiments on optically trapped cesium 
0| these conditions indeed turned out to be essential for 
the attainment of Bose-Einstein condensation. 

The molecules formed in recombination are likely to re- 
main trapped in our setup. A previous experiment shows 
evidence of trapping of the Cs2 dimer in a CC>2-laser 
trap |26| . Away from Feshbach resonances, the magnetic 
dipole moment of the dimer in its last bound state has 
the same value as the sum of the dipole moments of two 
single Cs atoms, so that the levitation trap can hold the 
molecules under the same conditions as the atoms. This 
offers an intriguing way to study atom-molecule interac- 
tions both in ultracold thermal clouds and Bose-Einstein 
condensates. Moreover, the Feshbach resonances allow to 
change the molecular magnetic moment in a controlled 
way and thus to separate the molecules from the atoms 
to prepare pure molecular samples. 
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